Detection of radio frequency magnetic fields using nonlinear magneto-optical rotation 
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We describe a room-temperature alkali-metal atomic magnetometer for detection of small, high 
frequency magnetic fields. The magnetometer operates by detecting optical rotation due to the pre- 
cession of an aligned ground state in the presence of a small oscillating magnetic field. The resonance 
frequency of the magnetometer can be adjusted to any desired value by tuning the bias magnetic 
field. We demonstrate a sensitivity of 100 pG/vTIz (RMS) in a 3.5 cm diameter, paraffin coated cell. 
Based on detection at the photon shot-noise limit, we project a sensitivity of 20 pG/VHz (RMS). 

PACS numbers: PACS. 07.55. Ge, 32.80.Bx, 42.65.-k 



I. INTRODUCTION 

Detection of small oscillating magnetic fields is the 
cornerstone of experimental techniques such as nuclear 
magnetic resonance (NMR), magnetic resonance imaging 
(MRI) , nuclear quadrupole resonance (NQR) Q and has 
been used in tests of physics beyond the standard model 
[3. Most atomic magnetometers (for example, see Refs. 

are designed to detect slowly varying magnetic 
fields and hence are not ideally suited for the aforemen- 
tioned applications. In recent work, Savukov et al |j| 
demonstrated a tunable, radio-frequency (RF) alkali va- 
por magnetometer, achieving a sensitivity of 20pG/VHz, 
with the sensor operating at 190°C. 

Here we present an RF atomic magnetometer based 
on nonlinear magneto-optical rotation (NMOR) arising 
due to the response of an aligned atomic ground state 
to a small RF magnetic field near the Zeeman reso- 
nance frequency. The Zeeman resonance frequency can 
be tuned to any desired value by adjusting the bias 
magnetic field, yielding sensitivity to signals of arbi- 
trary frequency. The measurement involves a single low- 
power light beam (~ 50 — 100 /iW), and based on pho- 
ton shot-noise limited polarimetry, achieves a sensitivity 
of about 20pG/\/Hz (RMS). The magnetometer has a 
bandwidth of about 50-100 Hz depending on light power. 
Based on experimentally observed signal-to-noise ratio, 
we demonstrate a sensitivity of lOOpG/VHz (RMS). De- 
spite somewhat lower sensitivity than reported in Ref. 
Q, for many applications, the magnetometer described 
here has the significant advantage that the sensor oper- 
ates much closer to room temperature (the highest tem- 
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perature used in this work was 48° C). Additionally, us- 
ing an aligned state rather than an oriented state pro- 
duces smaller external magnetic fields which can po- 
tentially have a back reaction on the sample of inter- 
est. Furthermore, the low power requirements and single 
beam arrangement facilitate the use of microfabrication 
techniques, promising for the development of compact 
portable atomic magnetometers 0. 

Such a magnetometer may find application in NQR as 
suggested in Ref. where the signal occurs at a fixed 
resonance frequency or in NMR spectroscopy where high 
spectral resolution is required to observe small splittings 
of NMR lines, due to, for example, scalar spin-spin (J) 
coupling between nuclei of the form JIi • I2. Such cou- 
plings c an y ield valuable information on molecular struc- 
ture [2| and can be difficult to access in high field 
environments where the absolute field homogeneity and 
differences in diamagnetic susceptibility limit the spec- 
tral resolution. Hence, recent attention has been given to 
performing such measurements in a low field environment 
using broadband, low transition-temperature supercon- 
ducting quantum interference devices (SQUIDs) [9| or 
inductive detection . As inductive detection becomes 
less efficient at low frequencies, the technique described 
in this letter offers the possibility of significant gains in 
signal-to-noise ratio without requiring cryogenics. 



II. RADIO FREQUENCY NMOR RESONANCE 

An alkali-metal vapor contained in a glass cell with 
anti-relaxation coated walls is placed in a z directed 
bias field Bo = Bqz,, corresponding to Larmor frequency 
CIl = 9^bBo where [Ib is the Bohr magneton and 
g w 2/(2/ + 1) is the Lande factor. Linearly polarized 
light propagating in the x direction with polarization in 
the z direction, tuned to the Dl [F = 2 ->■ F' = 1) 
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to zero and the magnetic field is 
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FIG. 1: Linearly polarized light, resonant with the Dl 
(F — 2 — *• F' = 1) transition, with polarization vector along 
Bo, produces an aligned ground state via optical pumping. 
Double headed vertical arrows indicate laser induced transi- 
tions between ground and excited states; dashed lines indicate 
transitions due to spontaneous decay. Ground state popula- 
tions, indicated by the solid black bars, are schematic only. 
A small RF magnetic field oscillating close to Ql, transverse 
to Bo, establishes coherences between neighboring Mf states. 
Inset: surface whose radius represents the probability for find- 
ing maximal projection of ground state angular momentum 
along a given direction (see, for example Ref. fl^l h 



transition, passes through the cell, optically pumping 
an aligned state, as illustrated in Fig. ^ We apply a 
small oscillating magnetic field transverse to the bias 
field, B x = B\ cos co r ft and we work in the regime where 
g^isBi <C ^ re i so that the RF magnetic field induces only 
ground state transitions with |AMp| = 1. 

We begin by assuming the light power is weak enough 
so that the saturation parameter relating the optical ex- 
citation rate to the ground state relaxation rate 



d 2 E 2 Vbeam 
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is small compared to unity. In Eq. QJ, d is the electric 
dipole matrix element of the optical transition, E is the 
light electric field, Tr> is the Doppler broadened width of 
the optical transition, 7 re ; is the ground state relaxation 
rate, V ce ii is the volume of the cell and Vbeam is the vol- 
ume contained within the intersection of beam and cell. 
When k <C 1 only the rank 2 (quadrupole) polarization 
moment is pumped by linearly polarized light. 

The oscillating RF magnetic field can be resolved into 
components co- and counter-rotating (parallel or antipar- 
allel to the direction of Larmor precession respectively) , 
each of magnitude Bi/2. Transforming to the co-rotating 
frame the counter-rotating component rapidly averages 
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In steady state, an equilibrium is reached between op- 
tical pumping of alignment along the z axis, precession 
around B' and relaxation, resulting in an aligned state 
tilted away from the z axis, as shown inset in Fig. ^ 
When uj r f = Ql, the z component in Eq. @ vanishes re- 
sulting in the maximum angle between the aligned state 
and the z axis. When we transform back into the lab 
frame, the tilted alignment processes about the z axis. 
The tilted alignment generates optical rotation through 
linear dichroism (see, for example, Ref. p^h maximal 
when the alignment is in the yz plane and none when 
it is in the xz plane, resulting in polarization rotation 
of the light beam that is modulated at to r f- The ampli- 
tude of the polarization rotation is linear in B\ in the 
range q[IbB\ -C jrei considered here, as we have verified 
experimentally. 

The description becomes slightly more complicated for 
conditions of high light power and light frequency de- 
tuned from optical resonance. Under these conditions, 
AC Stark shifts can lead to differential shifts of the 
ground state energy levels. In conjunction with preces- 
sion in the RF magnetic field, this results in alignment- 
to-orientation conversion (AOC) (see Ref. [14| and refer- 
ences therein) in the rotating frame and a splitting of the 
RF NMOR resonance as discussed briefly below. Doppler 
broadening can also lead to AOC effects, even for reso- 
nant light An additional high light power effect is 
the generation of hcxadecapole (rank 4) polarization mo- 
ments |l5j . We find experimentally that optimal sensi- 
tivity is achieved when the saturation parameter is close 
to unity, and density matrix calculations indicate that 
the hexadecapole contribution to the ground state po- 
larization is small compared to that of the quadrupole 
contribution for these conditions. 



III. EXPERIMENTAL SETUP 

A schematic of the experimental setup is shown in Fig. 
121 The measurements reported in this work were per- 
formed with an evacuated, paraffin-coated spherical cell 
(3.5 cm diameter) containing isotopically enriched 87 Rb 
(nuclear spin / = 3/2). The paraffin coating enables 
atomic ground-state polarization to survive tens of thou- 
sand wall collisions 16,|l7j, leading to ground-state po- 
larization lifetimes r = l/7 re z ~ 160 ms in a 10 cm diam- 
eter cell . The cell is placed inside a double- wall oven, 
temperature-controlled by flowing warm air through the 
space between the walls of the oven so that the optical 
path is unperturbed. A set of four nested /x-metal lay- 
ers provides a magnetically shielded environment, with 
a shielding factor of approximately 10 6 ^^|. Inside the 
innermost shield (cubic in profile) is a set of coils for the 
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FIG. 2: Schematic of the experimental setup. An evacuated, 
paraffin coated cell is placed inside a double wall oven. Tem- 
perature is controlled by flowing warm air through the space 
between the oven walls. A set of mu-metal layers provides 
a shielded environment and a set of coils maintains a stable, 
homogeneous magnetic field in the z direction. An additional 
coil generates a small oscillating magnetic field in the x direc- 
tion. Linearly polarized light from an external cavity diode 
laser passes through the cell and a balanced polarimeter mon- 
itors the polarization of the light as it exits the cell. 
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FIG. 3: a: Synchronously detected optical rotation and b: 
transmission spectra as a function of light frequency for a light 
power of 60 /jW and u r f = Ql- 



control of all three components of the magnetic field. Im- 
age currents in the magnetic shields create an "infinitely" 
long solenoid. The atoms traverse the cell many times 
during the course of one relaxation period, effectively av- 
eraging the magnetic field over the cell, leaving our mea- 
surements insensitive to field gradients |2C|. We apply a 
static magnetic field Bq in the z direction and a small 
oscillating magnetic field B\ cos(ui r ft) in the x direction 
(unless stated otherwise, B\ = 110 nG and Bo as 10 mG). 
A well collimated beam with diameter w 3 mm from an 
external-cavity diode laser, propagates in the x direction 
with polarization vector in the z direction. Unless stated 
otherwise, these measurements were performed with the 
light tuned to the center of the F = 2 — > F' — 1 tran- 
sition (henceforth referred to as optical resonance). The 
polarization of the light leaving the cell is monitored us- 
ing a balanced polarimeter and detected synchronously 
using a lock-in amplifier. Number density was deter- 
mined by monitoring the transmission of a low-power 
beam through the cell as a function of laser frequency. 
For the measurements reported here, the cell tempera- 
ture was 48°C and the measured number density was 
n = 7 x 10 10 (within 20% of that expected from the satu- 
rated vapor pressure at this temperature) , corresponding 
to approximately 1 absorption length for resonant light. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

In Fig. we plot the in-phase component of the syn- 
chronously detected optical rotation as a function of light 
frequency for co r f = Q^. For these data, the light power 
was 60 ^W(850/iW/cm 2 ). In Fig. 03 we plot the par- 



tially saturated transmission curve under the same exper- 
imental conditions. The background slope of the trans- 
mission curve is due to varying laser intensity as the 
diode laser feedback grating is swept. The largest op- 
tical rotation occurs for light tuned near the center of 
the F = 2 — > F' = 1 transition, similar to observations 
of non-linear Faraday rotation induced by a static mag- 
netic field [2l| • At the light powers for which we obtained 
optimal sensitivity on the F = 2 component, optical ro- 
tation on the F = 1 component was at least an order 
of magnitude smaller than that produced by the F = 2 
component. 

In the main panel of Fig. 01 we show the synchronously 
detected in-phase (stars) and quadrature (squares) com- 
ponents of optical rotation for light tuned to optical res- 
onance and incident light power of 40 fiW. Overlaying 
the in-phase (quadrature) component is a fit to a single 
absorptive (dispersive) Lorentzian. As mentioned previ- 
ously, under conditions of high light power and detuning 
far from optical resonance, differential AC Stark shifts 
can lead to a modification of the quantum beat frequency 
for different AMp = 1 transitions resulting in a splitting 
of the resonance, as shown inset in Fig. The over- 
all shift of the resonance compared to that shown in the 
main panel is because the bias magnetic field differed by a 
factor of 5. AC Stark shifts were discussed in some detail 
in Rcf. |f| in the context of NMOR with frequency mod- 
ulated light. It was found that magnetometric sensitivity 
was reduced when the resonance was split and hence we 
do not focus on this behavior any further. However, we 
point out that AC Stark shifts may provide some degree 
of optical tunability of the Zeeman resonance and we will 
explore this possibility in future work. 

In Fig.Et we plot Av, the half-width at half-maximum, 
of the in-phase component of the RF NMOR resonance, 
as a function of light power (the distance from the cen- 
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FIG. 4: In-phase (stars) and quadrature (squares), syn- 
chronously detected optical rotation signal as a function of 
the frequency of the applied oscillating magnetic field. Over- 
laying the data (solid lines) are fits to a single dispersive or 
absorptive Lorentzian. For these data, the light was tuned ap- 
proximately to the center of the F = 2 — > F' — 1 transition 
and the power was 40 /iW. Inset: Optical rotation signals for 
light power of 360 /iW, cell temperature of 20°C and optical 
detuning far to the red of the F = 2 — * F' = 1 transition. 
AC Stark shifts result in a splitting of the resonance. Units 
are the same as in the main panel. The overall shift of the 
resonance is different from that in the main panel because the 
bias field differed by a factor of 5. 



ter of the resonance to the extrema of the quadrature 
signal is also given by 5u). Overlaying the data is a 
linear fit with zero-power width Ai/o = 9.7 Hz. The 
intrinsic polarization relaxation rate 7 re ; is related to 
Ai^o via 'jrei = Itx^Vq |27| . Ground state relaxation in 
paraffin coated cells is typically dominated by electron 
randomization during collisions with the cell walls and 
through alkali-alkali spin exchange collisions (see for ex- 
ample [22l l2a ] and references therein). The relaxation 
rate for the latter process is given by |25j 



Ise ~ -zV SE v reX n 
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Here use ~ 2 • 10 ~ 14 cm 2 i s the spin-exchange cross sec- 
tion, and v re i = y // 8kT/TTfl is the average relative speed of 
the atoms, fi is the reduced mass. The factor of 1/2 in Eq. 
P|l represents the approximate nuclear "slowing down" 
factor appropriate for a spin 3/2 nucleus. For a density 
n = 7-10 10 cm~ 3 , Eq. (J3J gives7s£ = 27r-4.2Hz, roughly 
a factor of 2 smaller than the experimentally measured 
relaxation rate. We attribute the excess relaxation to 
collisions with the walls. 

In Fig. 0d we plot the amplitude of the RF NMOR 
resonance shown in Fig. 0] (defined as the maximum of 
the in-phase component) against the left axis as a func- 
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FIG. 5: a: Half-width at half-maximum of the in-phase com- 
ponent of the RF NMOR resonance shown in Fig. 0] as a 
function of light power. The solid line overlaying the data is 
a linear fit. b: The amplitude (squares) of the RF Zeeman 
resonance as a function of light power against the left axis. 
Against the right axis we plot the projected sensitivity (tri- 
angles) based on the amplitude of the resonance and photon 
shot-noise limited polarimetry. Points are joined as a guide 
for the eye. 



tion of input light power. The amplitude increases as a 
function of light power for low light power, until reaching 
a maximum at around 15/zW corresponding to k f=a 1.5. 
Beyond saturation the amplitude decreases due to light 
broadening. Against the right axis in Fig. [3Jd we plot the 
sensitivity of the magnetometer assuming a photon shot- 
noise limited polarimeter sensitivity S4> p h — V(2y^V0 
where is the number of photons per second exiting 
the cell. Light power is measured both before and af- 
ter the beam passes through the cell to accurately take 
into account optical absorption. Optimum sensitivity of 
about 20 pG/ %/Hz (RMS) is reached at about 40 - 50 //W 
input light power 28]. The bandwidth of the magne- 
tometer (the range of frequencies over which the signal is 
greater than half the value when u> r f — Ql) is given by 
the full-width at half-maximum, about 50 Hz at 40 fiW. 
By increasing light power to 100 /iW the bandwidth may 
be doubled with little loss in sensitivity. 

Atomic shot noise places a fundamental limit on the 
magnetometer sensitivity, given by 
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4pG/VHz(RMS) . (4) 



The factor of 2 is due to the fact that the magnetome- 
ter is sensitive only to the co-rotating component of an 
oscillating magnetic field. For an optimized magnetome- 
ter, one would expect that photon shot noise (adding in 
quadrature to the atom shot noise) would be compara- 
ble to the atomic shot noise [24| and hence Eq. (0J must 
be multiplied by a factor of y/2 for a fair comparison. 
Thus, for optimized light power, the sensitivity based on 
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photon shot noise in Fig. [5] is within a factor of 2 or 3 
of the fundamental limit. This indicates that the opti- 
cal pumping and probing scheme in this work is an ef- 
ficient method for detecting spin precession, despite the 
fact that atoms can be pumped into the optically dark 
F = 1 state. Some improvement in sensitivity may be 
achieved by increasing number density, however, in the 
regime where spin-exchange is the dominant relaxation 
process, r scales inversely with density so that the shot- 
noise limit is independent of density at sufficiently high 
densities. We point out that increasing the number den- 
sity can yield higher bandwidths. 

A complete theoretical treatment of NMOR is difficult 
because of the presence of hyperfine structure, Doppler 
broadening, velocity mixing and evolution in the dark. 
Following the method outlined in Ref. we performed 
a simplified steady-state density matrix calculation on an 
F = 2 — > F' = 1 transition which neglects these issues, 
but qualitatively reproduces the salient features of our 
experimental data. The hamiltonian is written in the 
rotating-wave approximation, neglecting terms counter 
rotating at either the optical or radio frequencies. The 
density matrix evolution equations are then formed, in- 
cluding terms describing spontaneous decay of the ex- 
cited state, and atoms entering and leaving the interac- 
tion region (transit relaxation), and solved numerically. 
For conditions of high light power and detuning far from 
optical resonance, the model reproduces the splitting of 
the RF NMOR resonance shown inset in Fig. 21 When 
the light is tuned to optical resonance, a single feature is 
observed in the RF dependence of the optical rotation. 
The calculated power dependence of the amplitude of the 
resonance is similar to the experimentally observed be- 
havior shown in Fig. [S] The model indicates that for 
light power that maximizes optical rotation (saturation) , 
the hexadecapole contribution to the ground state polar- 
ization is small, (roughly 10%) compared to that of the 
quadrupole contribution. 



V. RF MAGNETOMETER PERFORMANCE 

In Fig. we plot the noise spectrum as measured by 
an SRS770 spectrum analyzer at the output of the bal- 
anced polarimeter. The large peak is an applied field 
of 83 nG (RMS) to calibrate the magnetometer. Base- 
line noise is about 100 pG/vHz (RMS), falling some- 
what short of the sensitivity estimates based on pho- 
ton shot noise detection in Fig. [SJx Deformation of 
the beam by the cell resulted in a factor of roughly 
5 loss of light power, pointing towards an obvious im- 
provement for future work. In order to asses the per- 
formance of the polarimeter, shown inset in Fig. [f)| is 
the measured noise floor (squares) as a function of light 
power incident on the polarimeter. The dashed line rep- 
resents photon shot-noise 84> p h = l/(2 y JW^y l ) — 0.35 • 
/irad^/iW/ \/Hz (RMS). Polarimeter noise can be pa- 




6800 6900 7000 7100 7200 7300 
Frequency (Hz) 

FIG. 6: RF magnetic noise spectrum. The large peak is an 
applied field of 83nG(RMS) at Sl L . Light power was 40 /xW. 
Shown inset is the polarimeter noise (squares) as a function 
of light power incident on the polarimeter. The solid line is 
a fit base on Eq. Q and the dashed line represents photon 
shot noise. 



rameterized by 

H = sJc 2 ph /P + Q mp /P 2 - (5) 

Here P is the power incident on the polarimeter 
and C, p h and Camp parameterize photon shot noise 
and amplifier noise respectively. The solid line over- 
laying the data is a fit based on Eq. JSJ, result- 
ing in Camp = 0.55 ^rad/^W/VHz (RMS) and ( p h — 
0.41 jUrad^/juW/Hz (RMS), close to the theoretically pre- 
dicted value. Hence, amplifier noise is the dominant con- 
tribution for incident light power less than about 2 fiW 
and photon shot-noise dominates for higher light power. 



VI. CONCLUSION 

In conclusion, we have demonstrated a simple atomic 
magnetometric technique for the measurement of small 
RF magnetic fields based on a ground state Zeeman res- 
onance and detection of non-linear magneto optical ro- 
tation in an alkali-metal vapor. Based on photon shot 
noise detection we estimate a sensitivity of approximately 
20 pG/\/Hz (RMS) in a 3.5 cm diameter cell with a band- 
width of approximately 100 Hz. Estimates of the atom 
shot noise limit are within a factor of 2-3 of this limit, 
confirming that the present optical pumping and prob- 
ing scheme is an efficient method for probing spin pre- 
cession. Based on the actual signal to noise ratio, we 
have achieved a sensitivity of about 100 pG/ %/Hz (RMS). 
With several technical improvements, we anticipate a 
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factor of 3-5 gain in sensitivity. Optimization of num- 
ber density may yield some further gains in sensitivity 
as well as bandwidth. The magnetometer operates near 
room temperature, making it particularly attractive for 
applications in NMR. One possible such application is 
the measurement of a scalar, electron-mediated nuclear 



spin-spin coupling which can yield valuable information 
on molecular structure. The authors thank J. Higbie for 
useful comments and discussions. This work is supported 
by an ONR MURI program and KBN grant # 1 P03B 
102 30. 
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